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Abstract Resonant absorption of fast magnetoacoustic (FMA) waves in an inhomoge- 
neous, weakly dissipative, one-dimensional planar, strongly anisotropic and dispersive 
plasma is investigated. The magnetic configuration consists of an inhomogeneous mag- 
netic slab sandwiched between two regions of semi-infinite homogeneous magnetic 
plasmas. Laterally driven FMA waves penetrate the inhomogeneous slab interacting 
with the localised slow or Alfven waves present in the inhomogeneous layer and are 
partly reflected, dissipated and transmitted by this region. The presented research aims 
to find the coeflicient of wave energy absorption under solar chromospheric and coronal 
conditions. Numerical results are analyzed to find the coefficient of wave energy absorp- 
tion at both the slow and Alfven resonance positions. The mathematical derivations are 
based on the two simplifying assumptions that (i) nonlinearity is weak, and (ii) the 
thickness of the inhomogeneous layer is small in comparison to the wavelength of the 
wave, i.e. we employ the so-called long wavelength approximation. Slow resonance 
is found to be described by the nonlinear theory, while the dynamics at the Alfven 
resonance can be described within the linear framework. We introduce a new concept of 
coupled resonances, which occurs when two different resonances are in close proximity 
to each other, causing the incoming wave to act as though it has been influenced by 
the two resonances simultaneously. Our results show that the wave energy absorption 
is heavily dependent on the angle of the incident wave in combination with the incli- 
nation angle of the equilibrium magnetic field. In addition, it is found that FMA waves 
are very efficiently absorbed at the Alfven resonance under coronal conditions. Under 
chromospheric conditions the FMA waves are far less efficiently absorbed, despite an 
increase in efficiency due to the coupled resonances. 

Keywords: Magnetohydrodynamics (MHD); Sun; atmosphere; Sun; magnetic held; 
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1. Introduction 

The complicated interaction of the motion of plasma with the magnetic fields is one 
of the most interesting processes in the solar atmosphere. A highly non-uniform and 
dynamical system, such as the solar atmosphere, is a perfect medium for magnetohy- 
drodynamic (MHD) waves. These waves are able to transport momentum and energy 
along and across the magnetic field which can be dissipated, through a series of mech- 
anisms, one of the most efficient being resonant absorption. The existence of resonant 
absorption is due to the coupling of global oscillations to local waves in a non-ideal 
plasma. 

At resonance interacting dynamical systems can transfer energy to each other In 
this context, resonant absorption has been suggested as a mechanism to create supple- 
mentary heating of fusion plasmas, but was, however, later rejected due to technical 
difficulties (see, e.g. |Tataronis and Grossmann, 1973| [Tataronis and Grossmann, 1973[ 
IGrossmann and Tataronis, 1973', 'Grossmann and Tataronis, 1973", 'Chen and Hasegawa, 1974'| 
|Chen and Hasegawa, 1974 ; Hasegawa and Chen, 1976 Hasegawa and Chen, 1976 ). (lonson, 1978i )l 
suggested, for the first time, that resonant MHD waves may be a means to heat mag- 
netic loops in the solar corona. Since then, resonant absorption of MHD waves has 
become a popular and successful mechanism for providing some of the heating of 
the solar corona (see, e.;?. |Poedts eTa l., 1989l |Poedts et al., 1989) |Poedts et al., 1990a[ 
IPoedts et al., 1990a';'Poedts et al., 1990bi Poedtse Fal., 1990bl[Po edts et al. , 1990c| |Poedts et a l., 1990ct | 
[Ofman and Davila, 1995 Ofman and Davila, 1995]|Behen et al.,l999 , Belien et al., 1999j .| 

More recently resonant absorption has acquired a new applicability when the ob- 
served damping of waves and oscillations in coronal loops has been attributed to res- 
onant absorption. Hence, resonant absorption has become a fundamental constituent 
block of one of the newest branches of solar physics, called coronal seismology (which 
was first suggested by |Uchida et al., 1973[|Uchida et al., 1973t see also |Roberts et al., 1984H 
[Roberts et al., 1984] ). It has since been applied to explain the rapid damping of kink 
oscillations in coronal loops by numerous authors (see, e.g. [Nakariakov et al., 1999[ 
[Nakariakov et al ., 1999l |Aschwanden et al., 1999"Aschwanden et al, 1999VNakariakov and Ofman, 200T] | 
[Nakariakov and Ofman, 2 001 , R uderman and Roberts, 20 02 , R uderman and Roberts, 20021 1 
[Goossens et"aL72 002, Goosse ns et al., 2002[[Terradas et al.r2008|[Terradas et al., 2008t l 
[Terradas et al., 2010t |Terradas et al., 20T0| . In addition, resonant absorption has been 
proposed as a mechanism to explain the observed loss of energy of acoustic oscilla- 



tions in the vicinity of sunspots (see, e.g. Hollweg, 1988 HoUweg, 1988 [Lou, 1990| 


Lou, 19901 


Sakurai et al., 1991 Sakurai et al., 1991|[Goossens and Poedts, 1992"Goossens and Poedts, 1992[| 


Erdelyi anc 


Goossens, 1995 Erdelyi and Goossens, 1995; Goossens et al., 1995 Goossens et al., 1995^.| 



The efi'ect of resonant absorption has mostly been studied in linear MHD, mainly 
because the linear approach is easier to handle mathematically. However, linear theory 
shows that in the vicinity of the resonant position, the amplitude of oscillations can get 
large even when dissipative effects are taken into account resulting in a breakdown 



of the linear approach. The seminal study by ( [Ruderman et al., 1997d| l showed that 
this is exactly what happens at the slow resonance in isotropic plasmas. The ampli- 
tude of oscillations far away from the resonance are small, but close to the resonant 
position nonlinearity must be taken into account. The study was then extended by 
( [Ballai et al., 1998bl l to anisotropic plasmas characteristic for the solar corona. More 
recently, ( [Clack and Ballai, 2008| l expanded the theory to take into account dispersion 
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due to Hall currents. All analytical studies on nonlinear resonant absorption focussed 
on resonant slow waves despite Alfven waves being considered more important in the 
process of heating. Recently, (jCl ack et al., 2009| showed that the Alfven resonance 
can be described within linear MHD framework as the effects of largest nonlinear and 
dispersive terms identically cancel. 

Observations have shown that the solar atmosphere is not only an inhomogeneous 
and structured medium, but is also highly dynamic. Solar magnetic structures, such as 
coronal loops, sunspots and prominences, exhibit large-scale motions (flows). Modern 
high resolution satellites, however, are still not accurate enough to derive the detailed 
diagnostics of these flows (present day observations are limited to a resolution of ap- 
proximately 500 km, however, to fully resolve the processes creating these flows we 
could need a resolution of, conservatively, 1 km). On the other hand, theoretical studies 
have shown that these motions can considerably change the efficiency of MHD wave 



dissipation (see, e.g. Erdelyi and Goossens, 1996 Erdelyi and Goossens, 1996 1. The ef- 



fect of steady (shear) flows on the resonant behaviour of nonlinear slow waves was first 



discussed by (Ballai and Erdelyi, 19981 



Many studies of resonant absorption consider only the sound (or slow) and Alfven 
waves as excellent candidates for chromospheric and coronal heating (see, e.^.'Ruderman et al., 1997c'| 
Ruderman et al., 1997cl|BaUai et al., 1998a||Ballai et al., 1998a, Ruderman, 2000 Ruderman, 2000, | 



Vasquez, 2005 Vasquez, 2005[ l. However, Alfven waves can only carry energy along 



the magnetic field lines and slow waves are only capable of carrying 1 - 2% of their en- 
ergy under coronal (low plasma-j6) conditions. In contrast, fast magnetoacoustic (FMA) 
waves may have an important role in explaini ng the upper chrom ospheric and coro- 



nal temperatures as has been shown by, e.g. ( Cadez et al., 1997 1; ( |Csik et al., 1998| l 



and (Erdelyi and Ballai, 2001 1. FMA waves are magnetic waves which can propagate 
across the magnetic field lines. They are compressional so they are subject to dissi- 
pation by, e.g., viscosity, heat conduction. Landau and transit-time damping, in the 
high-frequency limit where the Coulomb collisions are ineffective. 

It is well known that in order to have acceptable heating by FMA waves coming 
from lower regions (e.g. generated by convective motions below the photosphere) the 
waves should not be reflected by the steep rise of the Alfven and / or slow wave speed 
with height, nor should they become evanescent. However, this is difficult to prevent. 
On the other hand, to interact resonantly with the slow waves, FMA waves would only 
have to propagate to the upper chromosphere, since the ratio between characteristic 
speeds (slow and Alfven) is much smaller than in the solar corona. Only the high- 
frequency waves of the full FMA spectrum can reach the corona (with periods of a 
few tens of seconds). The energy flux density of fast waves at the bottom of the corona 
required for significant heating is consistent with the upper limit on acoustic waves 



(Athay and White, 1978 Athay and White, 1978 1 



There is an alternative to the FMA waves propagating from the photosphere towards 
the solar corona. FMA waves (as well as slow and Alfven waves) may be generated 
locally in the chromosphere and corona by, e.g. phenomena involving magnetic re- 
connection. It was suggested by ( [Parker, 1988| l that shuffling the magnetic field in the 
solar atmosphere (by convective motion in lower regions) builds up magnetic stresses 
which can be released through, e.g. reconnection providing the energy to maintain the 
high temperature of the solar upper atmosphere. It has been shown that reconnection 
events can produce MHD waves (see, e.g. |Roussev et al., 200Ib|[Roussev et al., 200Ic[ 
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IRoussev et al., 2001a|[Roussev et al., 2001b|[Roussev et al., 2001c||Roussev et al., 200Ta| | 
through the process of nanoflaring. The required amount of nanoflares has not yet been 
observed due to their localised nature. Even extending the theory to larger observable 
scales (explosive events, micro-flares, blinkers, etc.) still does not prove this theory 
( |Aschwanden, 1999l|Aschwanden, 1999| l. Even though the energy released via recon- 
nection may not be adequate to provide the heating required, it is still assumed that 
the amount of magnetic energy built up by the shuffling can provide enough energy 
to heat the upper atmosphere. The necessary supplementary heating may come from 
reconnection-driven, locally generated MHD waves. 

The present paper has three major aims. First, is to study the nonlinear resonant 
interaction of FMA waves coupled into the slow continua in strongly anisotropic and 
dispersive plasmas (conditions typical of the solar upper chromosphere). Secondly, 
to numerically investigate the resonant absorption of FMA waves at the Alfven res- 
onance for conditions typically found in the solar corona. Finally, we numerically 
analyse the resonant absorption of FMA waves at a coupled (slow and Alfven) reso- 
nance in the solar chromosphere. The paper takes advantage of the analytical study by 
( |Clack and Ballai, 2009b| ) where they derived the coefficient of wave energy absorption 
for FMA, but did not discuss the implications of their findings in fuU based on concrete 
solar physics applications. 

The paper is organized as follows. The next section describes the equilibrium model 
and any mathematical and / or physical assumptions made. In Section 3 we give the 
definition of the coefficient of wave energy absorption at the slow and Alfven resonant 
positions. Section 4 displays our numerical results and our interpretation of the derived 
coefficients. Finally, in Sect. 5 we summarize, draw our conclusions and discuss the 
applicabihty of the study. 

2. The Equilibrium and Assumptions 

The magnetic field and density topology of the solar atmosphere is highly complex 
and, therefore, there are many degrees of freedom. However, to make our investiga- 
tion tractable mathematically, we study the interaction of incident FMA waves in a 
one-dimensional plasma. The dynamics and absorption of waves will be studied in 
a Cartesian coordinate system. The configuration consists of an inhomogeneous and 
dissipative magnetised plasma < x < jcq (Region II) sandwiched between two semi- 
infinite regions with homogeneous and ideal magnetised plasmas for x < Q and x > xq 
(Regions I and III, respectively). The equilibrium state is shown in Figure[T] which has 
been reproduced from ( |Clack and Ballai, 2009b| l. The equilibrium density and pressure 
are denoted by p and p. The equilibrium magnetic field, B, is unidirectional and lies 
in the yz-plane. In what follows the subscripts "e", "0" and "i" denote the equilibrium 
quantities in the three regions (Regions I, II, III, respectively). It is convenient to intro- 
duce the inclination angle, a, between the z-axis and the direction of the equilibrium 
magnetic field, so that the components of the equilibrium magnetic field are 

By - Bsina, B, - Bcosa. (1) 

All equilibrium quantities are continuous at the boundaries of Region II, and they 
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Figure 1. Schematic representation of the equihbrium state. Regions I (x < 0) and III (x > 0) contain a 
homogeneous magnetised plasma and Region II (0 < .v < xq) an inhomogeneous magnetised plasma. The 
shaded stiip shows the dissipative layer embracing the ideal resonant position Xc(a) ■ Image reproduced from 
Clack and Ballai (2009b). 



satisfy the equation of total pressure balance. It follows from the equation of total 
pressure balance that the density ratio between Regions I and III satisfy the relation 



Pe 

where the Alfven and sound speeds are 



24- 



VA = 




(2) 



(3) 



with y being the adiabatic exponent and jUo is the permeability of free space. 

We consider an equilibrium such that the plasma in Region III is both hotter and 
more rarefied than in Region I and assume a simple monotonic linear profile for all 
variables in Region II. This choice provides a single unique resonant surface for both the 
slow and Alfven resonances. The objective of the present paper is to study the resonant 
absorption of FMA waves at slow and Alfven dissipative layers. There are two cases 
which we wish to investigate. First, we study the absorption of fast waves at Alfven 
resonance as a possible scenario of the interaction of global fast waves (modelling FIT 
waves) and coronal loops. Secondly, we analyse the absorption of fast waves at the 
slow resonance as a possible scenario for heating in the upper chromosphere. When 
a fast wave interacts with a slow resonance, an Alfven resonance is also present. We 
consider this scenario in our analysis and study the total (i.e. coupled) absorption in the 
inhomogeneous layer. 

In an attempt to remove other effects from the analysis we consider the incoming 
fast wave to be entirely in the xz-plane, i.e. ky - 0. The resonant positions are located 
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where the global fast wave speed coincides with that of the local slow / Alfven speed. 
The dispersion relation for the impinging propagating fast waves takes the form 



1 



/t2 



{vl + cl)+ [vl + clf -Avici cos^ 



(4) 



where is the angle between the direction of propagation and the background magnetic 
field within the xz-plane and k = {k^ + k^y^^. For the sake of simplicity, the ratio k jk~ 
is denoted k^. so that the dispersion relation (|4| becomes 



k2 9 11 



vi + Cc 1 + 



1 



1/2 



(5) 



where 1 + kI - 1 / cos^ (/>. 

We assume the plasma is strongly magnetized in the three regions, such that the 
conditions Wi(e)Ti(e) » 1 are satisfied, where Wi(e) is the ion (electron) gyrofrequency 
and Ti(e) is the ion (electron) collision time. Due to the strong magnetic field, trans- 



port processes are derived from Braginskii's stress tensor (see, e.g. Braginskii, 1965 



Braginskii, 1965'; 'Erdelyi and Goossens, 1995' 'Erdelyi and Goossens, 1995', 'Ruderman et al., 1997aH 
Ruderman et al., 1997a , Ruderman et al., 1997b , Ruderman et al., 1997b; Mocanu et al., 2008t | 



Mocanu et al., 20(j8l l. As we deal with two separate waves (slow and Alfven), we will 



need to choose the particular dissipative process which is most eflicient for these waves. 
For slow waves, it is a good approximation to retain only the first term of Braginksii's 
expression for viscosity, namely compressional viscosity ( [Hollweg, 1985|[Hollweg, 1985| l.[ 
In addition, in the solar upper atmosphere slow waves are sensitive to thermal conduc- 
tion. In a strongly magnetized plasma, the thermal conductivity parallel to the magnetic 
field lines dwarfs the perpendicular component, hence the heat flux can be approximated 
by the parallel component only ( [Priest, 1984| [Priest, 1984| l. On the other hand, since 
Alfven waves are transversal and incompressible they are affected by the second and 
third components of Braginskii's stress tensor, called shear viscosity ( [Clack et al., 2009) 
[Clack et al., 2009| l. Finally, Alfven waves are efficiently damped by finite electrical 
conductivity, which becomes anisotropic under coronal conditions. The parallel and 
perpendicular components, however, only differ by a factor of 2, so we will only con- 
sider one of them without loss of generality. All other transport mechanisms can be 
neglected. 

The dynamics of nonlinear resonant MHD waves in anisotropic and dispersive plas- 
mas was studied by ( [Clack and Ballai, 2008| l and ( [Clack et al., 2009) 1. They derived the 
governing equations and connection formulae necessary to study resonant absorption in 
slow / Alfven dissipative layers. A further paper by ( [Clack and Ballai, 2009b[ ) studied 
the interaction of fast waves with the slow and Alfven resonances and derived the coeffi- 
cients of wave energy absorption. They intentionally overlooked the complication of an 
Alfven resonance being present at the same time as the slow resonance. In the present 
paper, however, we include both resonances. These conditions would be typical of the 
upper chromosphere where the plasma-yS goes from being larger than unity to becoming 
smaller than unity, retaining the monotonic dependence of characteristic speeds in the 
inhomogeneous layer. 
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The key assumptions we use in the present paper in order to facihtate analytical 
progress are connected to the strength of nonlinearity and the wavelength of the in- 
coming wave. First, from the very beginning we assume that nonlinearity, at the slow 
resonance, is weak. This simplification allows us to study the nonlinear governing equa- 
tion derived by ( |Clack and Ballai, 2008| l analytically. ( |Ruderman, 2000[ ) investigated 
the absorption of sound waves at the isotropic slow dissipative layer in the limit of 
strong nonlinearity. In their analysis nonlinearity dominated dissipation in the resonant 
layer which embraces the dissipative layer They concluded that nonlinearity decreases 
absorption in the long wavelength approximation, but increases it at intermediate values 
of kxo, however, the increase is never more than 20%. Secondly, we also assume that 
the thickness of the inhomogeneous region (Region II) is thin in comparison to the 
wavelength of the impinging fast wave, i.e. kxo 1, which has two implications. The 
first is to enable us to neglect terms of the order of k^x^ (and above) in the calculations. 
The second consequence is more subtle. To have a coupled resonance, the incoming 
FMA wave must resonantly interact with an Alfven resonance and be partially trans- 
mitted and then resonantly interact with a slow resonance. On one hand, if the Alfven 
and slow resonances are more than one wavelength, A: apart the transmitted wave will 
have decayed by the time it reaches the slow resonance. On the other hand, if the Alfven 
and slow resonance are within one wavelength of each other then the transmitted wave 
will be able to interact with the slow resonance. So, if an environment is set up such 
that both an Alfven and a slow resonance occur within the inhomogeneous region the 
distance between the resonances, xi , must satisfy the condition xi < xq. Hence, we have 
kxQ <«; 1 xo «: k'^ jc, «: k'^. 

We will use the same governing equations, absorption coefficients and results as 
obtained by ( |Clack and Ballai, 2009b) . We will only derive new results necessary to 
investigate the processes at the coupled resonance. 



3. Coefficient of Wave Energy Absorption 

The basic premise of our model is that one single monochromatic wave (the FMA wave) 
interacts with the inhomogeneous layer and one (or more) wave leaves it. The difference 
of the energy flux entering and leaving the resonance quantifies the energy resonantly 
absorbed inside the inhomogeneous region. This energy can then be converted into heat 
by dissipation. The coefficient of wave energy absorption is just the percentage of the 
energy entering the inhomogeneous region that is resonantly absorbed. 

( |Clack and Ballai, 2009b[ ) derived the coefficient of wave energy absorption of fast 
magnetoacoustic waves at both the Alfven and slow resonances, but the derivation 
requires a long and cumbersome mathematical treatment so we do not repeat their 
working, we just recall their main results. It should be noted that when they derived 
the coefficients of wave energy absorption they relied on connection formulae, which 



are assumed to be accurate ( Spruit and Bogdan, 1992 Spruit and Bogdan, 19921. At 



the Alfven resonance the wave dynamics are described by linear theory, as proved by 
( |Clack et al., 2009) 1. At resonance, the incoming monochromatic FMA wave is partly 
reflected, transmitted and dissipated. The reflected wave is monochromatic and the 
transmitted wave carries no energy because it is evanescent in Region III. The coef- 
ficient of wave energy absorption can be written as (see, e.g. [Ruderman et al., I997c[ 
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IRuderman et al., 1997c[ [Clack and Ballai, 2009b| [Clack and Ballai, 2009b| l 

where Ilin and Ilout are the normal components of the energy fluxes, averaged over 
a period, of the incoming and outgoing waves, respectively. Here is the plasma 
pressure in Region I and A„ denotes the total amplitude of the reflected wave at that 
order of approximation. Thus, the coefficient of wave energy absorption of FMA waves 
at the Alfven resonance is given by 



(Ta + lU.)^ + V- 



Ta = . , , + (7) 



where Ta, // and v represent 



POalAal ' pe(v2 - V^^COS^a)' 



(8) 



^-kV^ \ F-\x)dx. (9) 

A (y^ - I'Ai cos^ a\ Jo 

Here the subscript 'a' means an equilibrium quantity that has been calculated at the 
Alfven resonant surface, V is the phase velocity of the incoming FMA wave, given by 
Equation ([4]i, - + tan <p, Aa = -{dv^Jdx)^, cos- a and Ki is defined as 



\ ' (4 + vi,)(y2-4cos2a) 

The Cauchy principal part, is used in the definition of v because the integral is 
divergent. The function F{x) is the coefficient function of the highest order derivative 
in the Hain-Lust equation defining the continuous parts in the linear spectrum. It takes 
the form 

Po (v\ + cl](v- - cos^ a] (V- - cl. cos- a] 

F(x)^—^ -^-^^ '-^ (11) 

y4 - y2 (v^ + cl) + x'lcl cos2 a 

with cj - csVa(cj + v^)"'^^ being the cusp (tube) speed. 

The governing equation at the slow resonance have has been derived by ( [Clack and Ballai, 2008| l 
Similar to the Alfven resonance the incoming FMA wave is partly reflected, transmitted 
and dissipated. The transmitted wave is decaying in a similar way as before, however, 
the reflected wave includes higher harmonics in addition to the fundamental mode due 
to the nonlinearity at the resonance. Using the approximation of weak nonlinearity, in 
the first order of approximation, we find the fundamental mode of the reflected wave. 
In the second order approximation, we see that the higher harmonics are of the order 
of k^x^. In the third order of approximation, we find an additional contribution to the 
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reflected fundamental mode due to nonlinearity and dispersion. For even higher or- 
ders of approximation the contributions are no longer monochromatic, therefore, they 
are neglected (for full details on these estimations see, e.g. [Clack and Ballai, 2009b[ 
[Clack and Ballai, 2009b| ). The coefficient of wave energy absorption of FMA waves at 
the slow resonance is 

Tc - -5 T - C ; T + ^{i^ (12) 

where ^ and v take the same form as in Equations (jSj and (|9]l, and Tc is derived as 

nkV^ 



POcV^JAclcOS^Q' 



(13) 



Here the subscript 'c' means an equilibrium quantity has been calculated at the slow 
resonant surface and Ac = -{dc^ldx)aCos~ a. The first term in Equation ( [T2l i is the 
linear coefficient of wave energy absorption, and the second term is the nonlinear cor- 
rection (both terms are of the order of fcxo). The physical meaning of the small term ^ 
is buried within the nonlinear governing equation, but essentially describes the ratio of 
nonlinear and dissipative terms (for further details we refer to [Clack and Ballai, 2009b[ 
[Clack and Ballai, 2009b[ ). The small quantity ^ ensures that the contribution of non- 
linearity to the net coefficient of absorption is small. In addition, since the nonlinear 
correction is always positive, it is clear that the effect of nonlinearity and dispersion is 
to decrease the coefficient of wave energy absorption, a result already obtained earlier 
by ( [Ballai et al., 1998a| ). ( [Clack and Ballai, 2009bl i went on to show that when Hall dis- 
persion is also taken into account, the nonlinear term becomes 26 times larger, further 
reducing the overall coefficient. 

As mentioned earlier, the above discussed coefficients of wave energy absorption 
(taken separately) do not give the complete picture. In principle, every time a FMA 
wave interacts with a slow resonance an Alfven resonance is also present (it is easy 
to show that under coronal conditions, the vise-versa statement is not true). Usually, 
to avoid this happening we would align the magnetic field with the z-axis, and since 
d/dy = the Alfven resonance would vanish. However, to include Hall currents (which 
provides dispersive effects at the slow resonance) we must have an angle between the 
z-axis and the equilibrium magnetic field. In an attempt to tackle this problem, we 
investigate the interaction of FMA waves with an inhomogeneous region that contains 
both a slow and an Alfven resonance. To carry this out analytically, we must assume 
that the thickness of the inhomogeneous region is much smaller than the wavelength of 
the incoming wave, i.e. kxo <s: 1. We also assume that there is a single and unique slow 
and Alfven resonance (this assumption can be relaxed later). 

To find the coefficient of wave energy absorption at a coupled resonance the pro- 
cedure is identical to that carried out in ( [Clack and Ballai, 2009b| i, so we write out 
our findings rather than presenting all calculations in detail. First, we found that, at a 
coupled resonance, in the first order approximation (i.e. deriving the fundamental mode 
form), we cannot decouple the Alfven and slow reflected waves. This is in agreement 
with ( [Woodward and McKenzie, 1994a| l and ( [Woodward and McKenzie, 1994b| l who 
state that when Hall currents are present the slow and Alfven modes cannot always be 
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decoupled. Even though these studies were investigating stationary waves, we beheve 
the same applies here with resonant absorption, specifically, if we approximate that 
the wave interacts with both resonances simultaneously. The amplitude of the reflected 
wave in the first order of approximation is derived as 

A„ + Aci = -p j (14) 

T + 11 + IV 

where t — Ta +Tc, and yu and v are given by Equations ^ and Q, respectively. In higher 
order approximations we find that the Alfven resonance has no contribution to the 
amplitudes of the reflected waves. In particular, the addition to the fundamental mode 
due to nonlinearity and dispersion is identical to that found when no Alfven resonance 
is present, which is to be expected, since the dynamics at the Alfven resonance is linear 
and, therefore, should have no contribution to nonlinear effects. 

Using Equation (|6]l, the coefficient of wave energy absorption of FMA waves inside 
the inhomogeneous region is calculated to be (in the long wavelength approximation) 



(t + nY + 67r2 V^k^xl (fi^ + v^f 

In the next section, we will numerically analyse the coefficient of wave energy ab- 
sorption and, in particular, we investigate what is the significance of the angle of the 
incident wave and the inclination of the ambient magnetic field on the absorption of 
wave energy at the coupled resonance. 



4. Numerical Results 

Before proceeding to the mentioned analysis of the wave energy absorption, it is in- 
structive to consider the contributions to t due to the slow and Alfven resonances. 
We can easily calculate what percentage of t comes from each resonance by studying 
Equation ( fTSl l. Immediately it is clear that the nonlinear term comes wholly from the 
slow resonance, as we have already discussed. The linear term is a combination of the 
resonant absorption due to the slow and Alfven resonance. The percentage contributions 
of each resonance is found by calculating t^/t (% of slow contribution) and t-^/t (% of 
Alfven contribution). Using Equations ([8]), ( fT3] ) and ( fTSl l we obtain 

T K1+K2 T K1+K2 

where 

^1 = POalAalV"*, =POcVAclAc|sin^acos'a. (17) 

The particular form of Equation (fTSl l depends entirely on the choice of profile for 
the equilibrium quantities inside the inhomogeneous layer. For our analysis, in Region 
II (the inhomogeneous region) we have chosen a monotonically increasing linear pro- 
file for all equilibrium quantities (including characteristic speeds) of the form f{x) - 
/(O) + x/xo[f{xo) - /(O)]. Using the values vae = 28 km s"', cse = 34 km s"', v'Ai - 
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Figure 2. Left plot: percentage contribution to r in the coupled inhomogeneous layer due to the slow reso- 
nance in terms of the wave incident angle, (j), and the inclination angle of the magnetic field, a. Right plot: 
percentage contribution to t due to the Alfven resonance in terms if) and a. 

156 km s"', csi - 65 km s"' and pe = 5 x 10"" kg m""' for characteristic speeds 
we can obtain the variations given by Figure [2] The left plot in Figure |2] shows the 
percentage contribution to the resonant absorption due to the slow resonance, while 
the right plot shows the same for the Alfven resonance. It is easily seen that (f> does 
not have much effect on the contributions, however, a has a marked effect. At ff = 0, 
there is no contribution from the Alfven resonance - which is to be expected since, for 
this particular inclination, an Alfven resonance does not exist. As the angle a increases 
the contribution due to the Alfven resonance increases rapidly (and accordingly the 
contribution from the slow resonance decreases), and hy a - njA the Alfven resonance 
contributes 80% of the absorbed energy. At a = njl the slow resonance disappears, and 
so its contribution drops to zero. 

4. 1 . Alfven Resonance: Modelling the Interaction of BIT Waves with Coronal Loops 

Fast waves that are generated by the convection motion and propagate upwards are 
reflected by the strong density gradients in the upper chromosphere, so just a tiny 
proportion of FMA waves are able to reach the corona. Fast waves, however, can be 
generated in the corona by, e.g. flaring processes, in particular coronal mass ejections 
(CMEs). Global disturbances generated by CMEs, known as BIT (Bxtreme ultravi- 
olet Imaging Telescope) waves, are believed to be FMA waves ( |Ballai et al., 2005] 
[Ballai et al., 2005| ) propagating in the quiet Sun. In their propagation, BIT waves in- 
teract with active region loops setting them into motion (see, e.g. [Ballai et al., 2008| 
[Ballai et al., 2008| ). The present section is devoted to the study of resonant coupling 
of FMA waves (modelling BIT waves) with local Alfven waves in coronal loops. One 
major problem identified with the model we would like to use is that the BIT waves 
have very low speeds (300 - 500 km s"') compared with the accepted Alfven speed 
outside coronal loops (> 1000 km s"'). The problem associated with the speed of BIT 
waves could be resolved by assuming a steady rise in temperature and magnetic field 
strength as an active region is approached. The rise in temperature increases the sound 
speed, while the increase in strength of magnetic field will cause the Alfven speed to 
grow. Hence, the FMA wave speed will increase as it is a combination of the sound and 
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Alfven speeds. As the gradients are not steep enough to cause shocks or reflection, the 
EIT waves are accelerated to the local speeds as they approach active regions. 

Another assumption our model relies on is that vai > vab, ' e. the Alfven speed in Re- 
gion I is less than that in Region III. Since we can only integrate F{x) when we assume 
monotonic functions inside the inhomogeneous region we have to chose whether these 
functions will increase or decrease. If the functions decrease in the inhomogeneous 
region there will be no absorption since the FMA waves will never resonantly interact 
with the local Alfven waves. However, if the functions increase, we do achieve resonant 
interactions of the two modes. A more accurate model of reality at a coronal loop's 
edge would be a complicated non-monotonic function that both increases and decreases 
the local equilibrium parameters which could create resonant interactions. Essentially, 
there would be several resonant positions inside the inhomogeneous region, however, 
we cannot analytically solve such a model at present (due to the form that F(x) would 
take). Therefore, to provide a valuable insight into the efficiency of resonant absorption 
at coronal loops, we choose an unrealistic but tractable model. Mathematically, the 
above restriction can be easily represented by the condition that max[vA(x)] > va^ 
where x € (0, xq]. In addition, for a single unique resonance we need to impose the 
constraint dv\{x)/dx i^Q, x e (0, xq). It is also clear, from the outset, that under coronal 
conditions (y6 <k 1) FMA waves will never resonantly interact with local slow waves. 

To simulate conditions typical of the solar corona we consider VAe = 1200 km s"', 
cse = 200 km s"', = 1400 km s"', csi = 250 km s"' andpe = 1.33 x 10"'^ j^g m"^ 
[pi being calculated automatically via Equation (|2]i]. We select k - 5 x 10^^ m"' such 
that the incoming FMA wave has a period of about 100 s (consistent with the order 
of magnitude of the period of observed EIT waves in the solar corona). It should be 
mentioned that the analysis presented in Figure [3] is valid for any k as long as the 
dimensionless quantity Lyq satisfies the condition kxo <s: 1 . 

Figure|3]gives the wave energy absorption coefficient for kxQ - 0.01, 0.1, 0.5 and 1.0 
(the last one not being full consistent with our assumptions, but illustrates the trend). 
The angle of the incoming FMA wave, 0, takes values between and n/2, while the 
angle of the equilibrium magnetic field, a, only varies between and 7r/4, because be- 
yond this point the integrals calculated are divergent and the numerical analysis cannot 
resolve the Cauchy principal part. The coefficient of wave energy absorption should 
take a value between and 1 , and can be thought of the percentage of incoming wave 
energy transferred to local Alfven waves by resonance. We recognize that the plasma 
is unstable when the value of the coeflicient of wave energy absorption is negative, 
provided a flow exists. The plasma can also create over-reflection, and we observe this 
in the coefficient of wave energy absorption when its value becomes greater than unity. 
Both the phenomena of over-reflection and instabilities occur for a > n/4. 

The first thing to notice, about Figure[3j is that as kxo increases so does the coefficient 
of wave energy absorption. In the same way, the area (in the a, (p plane) over which 
absorption can take place is also increasing. In general, we can state that the most 
efficient absorption occurs at angles of inclination of which are larger than n/A, while 
the most efficient absorption with respect to a changes depending on kxQ. The larger 
the value of kxo the less influential a becomes on the most efficient absorption regions. 
At approximately a - tt/S a; 0.63 the absorption drops to zero for all values of kxo and 
(p, which could be explained by the inhomogeneous layer becoming transparent to the 
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Figure 3. The wave energy absorption coefficient of FMA waves at the Alfven resonance. Here we have 
VAc = 1200 kms"', csc = 200 km s"', va, = 1400 km s"', cs, = 250 km s"' and pe = 1-33 X 10"'^ kg m"^ 
The dimensionless quantity, kxo, takes the values (from top left to bottom right) 0.01, 0.1, 0.5 and 1.0, 
respectively. 

incoming wave, a phenomenon that cannot be explained by the present model (further 
study is needed to find out whether it is a numerical artefact or a physical property). 

It is also clear from Figure[3]that FMA waves are absorbed efficiently at the Alfven 
resonance, which is encouraging when thinking about BIT waves within the solar corona.| 
The BIT waves could be absorbed by Alfven resonances present in / or near coronal 
loops (arcades). In some cases all of the energy of the incoming wave can be absorbed, 
and dissipated by, e.g. viscosity. The variation in absorption due to combinations of a 
and (p may help explain why when an incoming wave impacts a coronal arcade, some 
loops oscillate more than others and some loops dim and while others get brighter 
When the frequency of the incoming FMA wave does not fall within the frequencies 
of the Alfven continuum the energy of the incident wave is, likely to be, transferred 
to the coronal loop as kinetic energy, thereby setting the loop into oscillation. These 
oscillations are studied in the framework of coronal seismology for the purposes of 
field and plasma diagnostics. 

We note here that, if the density is varied (within reasonable parameters), the absorp- 
tion rate is changed only slightly. If we change the density so it is really high or low 
(for the solar corona) the absorption starts dramatically changing, eventually leading 
to a breakdown of the numerical analysis. We do not show any plots of the variation 
of the absorption coefficient with density because the values at which the coefficient 
of wave energy absorption is noticeably changed are not consistent with observation 
of the solar corona (or chromosphere). The last significant variables to discuss are the 
equilibrium wave speeds. The characteristic wave speeds used here have been selected 
to be consistent with the environment of the solar corona. However, we know that there 
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is a plethora of wave speeds available in the inhomogeneous plasma of the corona. We 
have tested different values of equilibrium speeds, and the overall pattern of absorption 
is identical to that found in Figure |3] and the absorption coefficient remained relatively 
similar to those discussed above. 

4.2. Coupled Resonance: Modelling Chromospheric Absorption 

To model a coupled resonance, we consider a slow and an Alfven resonance so close 
together that the local waves can both interact with the incoming FMA wave. This 
means that we will lower our applicability region to the denser chromosphere. To match 
conditions typical of the chromosphere we adopt the values VAe = 28 km s cse - 
34kms"', VAi = 156 km s"', csi = 65kms"' and pe = 3.99 x 10"" kg m^l The 
specific values of equilibrium quantities can be changed, however, the overall variation 
trend of the coefficient of absorption will remain the same as shown here. The nonlinear 
correction to linear absorption at the coupled resonance comes solely from the slow 
resonance (as the Alfven resonance can be described by linear theory) and is very small 
when compared to the linear absorption coefficient (in line with the weak nonlinear limit 
imposed in the derivations). In Figure|4]we display the linear coefficient of wave energy 
absorption for the slow (top left panel), Alfven (top right panel) and coupled resonance 
(bottom left panel) for comparison together with the nonlinear correction to the coupled 
resonance (bottom right panel) for values of kxo = 0.5. We have chosen to only display 
one value of kxo because the plots for other values of kxo look similar to this one, the 
only difference is the magnitude of the coefficient of wave energy absorption where, in 
general, the greater kxo the higher rate of absorption. 

We have truncated the abscissa at a = 67r/25, because the integrals become divergent 
and the numerical analysis cannot resolve the Cauchy principal part. In Figure |4] the 
nonlinear coefficient of wave energy absorption is incredibly small since the ratio of 
magnetic and plasma pressures (the plasma-/?) has a value smaller than unity. When 
yS «: 1 the plasma pressure is low and, by extension, [examine Equation (fTST ll so is 
the nonlinear coefficient of wave energy absorption. In addition, the nonlinear coef- 
ficient must be multiplied by the very small parameter Therefore, the nonlinear 
correction of wave energy absorption at coupled resonances in the chromosphere is 
truly tiny in comparison to the linear wave energy absorption and acts to decrease the 
total coefficient of wave energy absorption. 

In general, we can state that the greatest linear wave energy absorption at the slow 
resonance occurs at small values of (0 < < 0.3) and small to moderate values of a 
(0.1 < a < 0.6), whereas the the greatest linear wave energy absorption at the Alfven 
resonance occurs at moderate values of a (0.6 < a < 0.8) and a wide range of values of 
<p (0.2 < (p < 1.2). For all values of kxQ a larger percentage of energy is absorbed at the 
Alfven resonance compared to the slow resonance, which is indicated by the coefficient 
of wave energy absorption having greater values at the Alfven resonance. We can also 
state that increasing kxo produces more absorption (of course up to the point where the 
mathematical model becomes invalid). 

The most interesting feature of Figure|4]is the coupled resonance. When there is no 
absorption at both the slow and Alfven resonance, the absorption at the coupled reso- 
nance is also zero and if there is absorption at only one of the slow or Alfven resonance 
the coefficient of wave energy absorption at the coupled resonance is identical to the 
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Figure 4. Comparison of tlie linear absorption at tlie slow (top left), Alfven (top right) and coupled (bottom 
left) resonances for kxo = 0.5. We also show the nonlinear absorption for the coupled (and slow) resonance. 



value at the single resonance. However, more interestingly, when there is absorption 
at both the single Alfven and slow resonance the absorption at the coupled resonance 
is always greater than the sum of the two single absorption coefficients. The greater 
coefficient of wave energy absorption at the coupled resonance compared to the two 
single resonances implies that there is more energy available for heating the plasma. 
We shall use Figure |4] to clarify what we have just explained. At the values of, e.g. 
a - 0,0 = 7r/2 in all three plots we can clearly see the coefficient of wave energy 
absorption is zero. When, e.g. a -Q,(p - 0.4 the coefficient of wave energy absorption 
at the Alfven resonance is zero, while at the slow and coupled resonance it is 0.11. 
If we change the angles such that, e.g. a - 6tiI 25, (p - 1.0 the coefficient of wave 
energy absorption at the slow, Alfven and coupled resonances are 0.067, 0.13 and 0.23, 
respectively, hence it is clear that 0.23 > 0.13 + 0.067 = 0.197. The same procedure 
can be carried out at any values of a and (p where there is absorption at both the single 
Alfven and slow resonances. 

There is also clear evidence in these plots that the efficiency at which FMA waves 
are absorbed at the Alfven resonance under chromospheric conditions is far lower that 
in the coronal counterpart. 



5. Conclusions 

In the present paper, we have numerically investigated the absorption of fast magnetoa- 
coustic (FMA) waves at individual and coupled slow and Alfven resonances, using the 
theory developed by ( |Clack and Ballai, 2009b[ ) where the coefficients of wave energy 
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absorption were derived analytically by applying the long wavelength approximation 

(kxQ <«; 1). 

We have shown that the absorption of incoming FMA waves depends heavily on 
the combination of the angle of incidence of the wave (0) and the angle of incUnation 
of the equilibrium magnetic field (a). At some combinations of (p and a, the quantity 

< which reduces absorption to zero, because the incident FMA wave can leak 
passed the resonances. Normally, the inhomogeneous layer is translucent to waves, 
but occasionally when the conditions are right, the inhomogeneous layer can become 
transparent so that they can pass through without undergoing resonant absorption. 

We introduced the concept of coupled resonances which could be important for the 
heating of the solar upper atmosphere, because FMA waves propagate throughout the 
solar atmosphere. When FMA waves interact with local waves at a slow resonance, an 
Alfven resonance is always present as well (as long as there is an angle between the 
magnetic field and the direction of wave propagation). We believe this is the most likely 
form of resonant absorption in the upper chromosphere. In the corona, the Alfven (and 
hence FMA) waves speeds are much larger then the slow speed, therefore, a laterally 
driven incoming FMA wave can encounter an Alfven resonance only. The governing 
equations for the slow and Alfven resonances are derived considering nonlinearity, and 
even though nonlinearity only slightly decreases the resonant absorption, in comparison 
to linear theory, it provides a device by which the heated plasma can be transported; 
the mean shear flow, creating turbulence, which can distort the inhomogeneous layer 
(see, e.g. 'Ofman and Pavila, 1995| [Clack and Ballai, 2009a| |Ofman and Davila, 19951 
[Clack and Ballai, 2009^ , enhancing absorption and transporting heated plasma away 
from the resonance. 

We assumed, for simplicity, that the equilibrium quantities inside the inhomogeneous 
layer increased monotonically, which, in reality, is not always the case. If we allow the 
equilibrium quantities to vary non-monotonically, the situation is changed. If we still 
use the long wavelength approximation, the calculations are almost identical to the ones 
produced in the present paper, however, instead of one position for the Alfven and slow 
resonances (different for each one) there could be several. Again, we would not be able 
to separate the individual outgoing waves, but the collective outgoing wave would be 
defined exactly as Equation ([T4l i. where r would equal the addition of all the t^s and 
TqS associated with the different resonant points. The form of the r^s and t^s would 
be identical to Equations ([8]l and fT3i . respectively, though the quantities would take 
diff'erent values corresponding to the different resonant positions. In general, this would 
produce greater absorption of wave energy inside the inhomogeneous layer, creating 
further heating possibility. 

Extreme ultraviolet (EUV) observations of active regions showed that, when coronal 
arcades oscillate under the influence of an external driver, some of the loops oscillate 
more than others, some of the loops do not oscillate at all and some become dimmer in 
the wavelength that they have been observed in. Part of this behaviour can be explained 
by the varying strength of the magnetic field inside the loops, however, this is not fully 
satisfactory. We speculate that some of this phenomenon can be explained by resonant 
absorption. The varying conditions allow for varying degrees of resonant absorption. 
The loops that oscillate the most do not have the conditions necessary for resonant 
absorption, so the incident waves just transfer their kinetic energy to the loops. Other 
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loops have the right conditions for resonant absorption, and some even have the condi- 
tions for multiple resonant positions. These loops will oscillate less, because the wave 
energy is being absorbed. The loops with the conditions for multiple resonances will 
oscillate the least (possibly not at all) as more and more energy of the wave is deposited. 
These loops should either get brighter or dimmer, since the loops are only observed in 
a single wavelength. If a loop brightens, the plasma is emitting more intensely in the 
filters wavelength, but if a loops dims the loop could be getting either cooler or hotter, 
as the plasmas emissions move out of the filters range. This explanation covers all the 
observed properties of some coronal arcade oscillations, and could explain why they 
are damped so quickly. The speculation must be viewed with caution as the heating due 
to absorption is likely to take place over a small area, in comparison to the width of the 
loop. We also add, however, that much more investigation is needed to substantiate this 
idea. 
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